Direct Measurement of the Free Energy of Aging Hard-Sphere Colloidal Glasses 



Rojman Zargar, 1 Bernard Nienhuis, 2 Peter Schall, 1 and Daniel Bonn 1 

1 Van der Waals-Zeeman Institute, Institute of Physics (loP) of the Faculty of 
Science (FNWI) University of Amsterdam, 1098 XH Amsterdam, the Netherlands 

2 Institute for Theoretical Physics, Institute of Physics (IoP) of the Faculty of 
Science (FNWI) University of Amsterdam, 1098 XH Amsterdam, the Netherlands 

The nature of the glass transition is one of the most important unsolved problems in condensed 
matter physics. The difference between glasses and liquids is believed to be caused by very large 
free energy barriers for particle rearrangements; however so far it has not been possible to confirm 
this experimentally. We provide the first quantitative determination of the free energy for an aging 
hard-sphere colloidal glass. The determination of the free energy allows for a number of new insights 
in the glass transition, notably the quantification of the strong spatial and temporal heterogeneity in 
the free energy. A study of the local minima of the free energy reveals that the observed variations 
are directly related to the rearrangements of the particles. Our main finding is that the probability 
of particle rearrangements shows a power law dependence on the free energy changes associated 
with the rearrangements, similarly to the Gutenberg-Richter law in seismology. 
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The simplest system that has become a paradigm for 
molecular glasses is a colloidal hard-sphere system in 
which spheres interact only through an infinite repulsion 
on contact pQ. The phase behavior of such hard-sphere 
systems is uniquely governed by the volume fraction <f> 
[2]: There is a liquid phase for <j> < 0.494, a crystalline 
phase for <f> > 0.545 and coexistence between the two 
for 0.494 < <f> < 0.545. If crystallization is avoided, 
for 0.545 < 4> < 0.58 the system is a supercooled liq- 
uid whereas for 4> > 0.58 the dynamics of the system 
becomes so slow that it is a glass [3HZ]- 
Because of the absence of interaction energies, the phase 
behavior of hard-sphere systems is dictated only by en- 
tropic contributions, implying that the statistical geome- 
try provides a direct route to study the thermodynamics 
of the system [8]. The usual way to do this is to de- 
termine the available volume, i.e. the volume available 
to insert an additional sphere into the system; there are 
exact relations that relate the available volume to the 
thermodynamics of the system [TU] . This method was 
successfully applied to colloidal liquids, but is unfortu- 
nately not applicable to glasses, because the cavities be- 
come extremely small at high densities and so the deter- 
mination of the available volume becomes prohibitively 
difficult [TU]. 

An alternative method relies on the use of the related free 
volume [llj , the volume over which the center of sphere 
can move when all other particles are held fixed. Unlike 
the available volume, the free volume is a local quantity 
that can be obtained for each individual spherical parti- 
cle. Since the free volume is typically much larger, it can 
be measured with much greater accuracy. 
In this Letter, we apply the free volume method to mea- 
sure the free energy of a colloidal glass. In this approach 
the system is divided into N cells constructed around 
the particle positions; choosing the cells so that they are 



singly occupied, the free energy of a hard-sphere system 
for a given configuration of the particles can be expressed 
directly in terms of only the free volume as [13] : 

N 
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in which Vf i is the free volume for particle V in its cell, 
and A = (/i 2 /27rmfcsT) 1 / 2 is the thermal wavelength; this 
form supposes there are no interactions other than hard- 
sphere. 

We determine Vf experimentally in a colloidal hard- 
sphere system; we use sterically stabilized fluorescent 
poly-methylmethacrylate (PMMA) particles suspended 
in a refractive index and density matched mixture of cis- 
decalin and cycloheptyl bromide (CHB) which are the 
best model system for hard-spheres [J2]. Our particles 
have a diameter a — 1.5fim and polydispersity ~ 7% to 
prevent crystallization. The crystal is obtained in a sec- 
ond system with a = 1.7 fim and polydispersity ~ 4%. 
In both systems an organic salt, TBAB (tetrabutylam- 
moniumbromide) , is used to screen any possible residual 
charges. Suspensions at volume fractions between 0.54 
and 0.62 are prepared by diluting sediments that are cen- 
trifuged to random close packing (<j) rcp ~ 0.64); for a dis- 
cussion on the accuracy of this procedure, see Ref. [13] . 
In order to verify the hard-sphere nature of the colloids 
the crystallization density, a very sensitive measure for 
deviations from the hard-sphere behavior, is measured: 
The results agree to that of the true hard-spheres within 
a fraction of a percent. To initialize the glassy samples 
we insert a small stirring bar into the chamber of each 
sample to stir and have a reproducible initial state for the 
sample at the start of the experiments. The initial time 
t = is defined as the beginning of the measurements 
that we start 20 minutes after ending the stirring. 
Using a fast confocal microscope (Zeiss LSM live) we ac- 
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FIG. 1. Distribution of the Voronoi volumes for the crystal 
and the glass at <j> = 0.60. 



quire 3D images of the system in a 105/ito x 105/xm x 
80/im volume to find the positions of ~ 2.5 x 10 5 individ- 
ual particles in three dimensions (Fig. [2](d)) . We do so for 
both disordered (liquid, glass) and ordered (crystalline) 
systems. Using the positions of particles, a Voronoi tes- 
sellation is generated [15] : A Voronoi tessellation divides 
the space into distinct, non-overlapping convex polyhe- 
dra which surround each point. The Voronoi volume is 
given by the polyhedron around the center of a sphere V 
consisting of the region of space closer to V than the cen- 
ter of all the other spheres. We find that the distribution 
of the Voronoi volumes for the particles in the crystal 
is much sharper than in the disordered system (Fig. [T]) 
indicating the same Voronoi volume for the particles in 
the crystal. In addition, the measured distribution for 
the Voronoi volumes agrees very well with simulations 

dang. 

The free volume however, is defined as the volume over 
which the center of sphere can move when all other par- 
ticles are fixed. Using the position of particles, we can 
then calculate the free volume for each particle which is 
assumed to be confined to its Voronoi cell. Following Ref. 
[13] and consistent with Ref. [TS] [TJ5] , the free volume as- 
sociated with a Voronoi cell is defined as the volume of 
a smaller cell generated from the Voronoi cell by moving 
the faces normally inside of a distance cr/2, where a is 
the particle diameter. It should be noted that the as- 
sumption of considering that the particles are confined 
in cells is valid only at high densities where the particles 
are so densely packed that no rearrangements happen on 
the time scale of obtaining the configuration. 
The observed distribution for the free volumes in the 
crystals is consistent with a Gaussian distribution (Fig. 
[2|a)). For the glass however, the observed distribution 
for the free volumes shows an exponential tail (Fig. [2^b) ) 
in line with random nature of the glassy state. For the 
same volume fraction, the average free volume for parti- 
cles in the glass is much smaller than that in the crystal 
(Fig. [5Ja),(b)), leading to a larger free energy for the 
glass; this is of course expected, because the crystal is 
the thermodynamically stable phase. 
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FIG. 2. (a),(b) Distribution of the free volumes at cj> — 0.54 
and (j> = 0.60. (a) Crystal; solid curves show Gaussian fit. (b) 
Disordered system; solid curves show exponential fit. (c) The 
free energy density F/V in units of k B T/nm 3 as a function of 
the volume fraction. Circles denote our free energy results for 
the crystal at volume fractions 0.54, 0.56, 0.58, 0.60 and 0.62 
that agree with the Hall [2UJ free energy (dashed curve). Stars 
show the free energy measured for the disordered systems at 
volume fractions 0.54, 0.56, 0.58, 0.60 and 0.62 at the begin- 
ning of the measurements while squares indicate the measured 
free energy after 5 days. The results agree rather well with 
the standard Carnahan-Starling [22] free energy (solid curve) . 
Crosses are the measurements of Dullens et al. |10j which are 
also compared to the Carnahan-Starling result, (d) Typical 
images of the confocal microscope and a part of the 3D recon- 
struction of the particle positions. Left: Crystal at <j> = 0.54, 
Right: Glass at <f> = 0.60. 



The colloidal crystal is used as a reference system for 
our measurements since there are other very accurate 
methods to determine the free energy for ordered systems 
[20l |2"T] . We find that the measured free energies for the 
hard-sphere crystal are in excellent agreement with the 
Hall equation of state (dashed curve in Fig. j2f;). On the 
other hand, the free energy for the disordered system for 
volume fractions up to the glass transition are found to be 
in good agreement with the standard Carnahan-Starling 
equation of state which is known to be very accurate for 
the liquid phase of the hard-sphere systems [22 . The 
free energy results of Dullens et al. [10 , who were only 
able to measure the free energy for <f> < 0.40 because they 
were using the available volume method, also agree with 
the Carnahan-Starling results. 

This method then allows, for the first time, to inves- 
tigate the free energy and its variations as function of 
both space and time for a glass. In Fig. [3^a) we show 
that, while the free energy of the crystal is constant as 
expected for an equilibrium system, it decreases signifi- 
cantly over a long time in disordered systems, signaling 
aging of the system [53] [M] • This is what is commonly 
assumed for the time evolution of non-equilibrium glassy 



3 




Time (min) 



FIG. 3. (a) Time evolution of the free energy. Free energy 
per particle is plotted versus time: Crystal (lower part) and 
disordered system (upper part) at different volume factions. 
'D' stands for disordered and 'C for crystal, (b) Aging of 
disordered systems. Free energy per particle is plotted versus 
time for <f> — 0.54 (brown) and (f> — 0.58 (black). The noise 
shows the effect of the local minima of the free energy. Note 
that the data in (a) corresponds to the very beginning of these 
curves, (c) Heterogeneity in the free energy distribution in a 
glass at (j> = 0.60 at time t — 82min . The color code shows 
the free energy per particle. The background indicates the 
free energy calculated for the whole system, (d) The free 
energy per particle versus time in one of the boxes considered 
in (c). 



systems, but we are the first to observe the temporal free 
energy decrease experimentally. For hard-sphere colloids 
it is well known that aging can be observed on similar 
time scales in e.g. the correlation function of the particle 
displacements [25| [26] . 

To check the effect of the particles size distribution on the 
calculated free energy, we add size polydispersity into the 
crystal as well as the glass artificially: A polydispersity 
of 7% distributed randomly in the glass decreases the 
free energy by 6% while a polydispersity of 4% randomly 
distributed in the crystal increases the free energy by 
an amount of 4%. This is much smaller than the typical 
fluctuations in the free energy and the difference between 
ordered and disordered systems that we consider here. 
To see whether there is a clear signature of glassy behav- 
ior, we measure the free energy for several days for a dis- 
ordered system with (j> = 0.54 and for </> = 0.58, which are 
commonly believed to be a supercooled liquid and a glass 
(Fig. [3^b)) [3J. In the supercooled liquid the free energy 
saturates to a horizontal plateau. For the glassy system 
at (j) — 0.58, after the quench into the glassy phase, the 
free energy evolves very slowly and indeed the glass never 
equilibrates on the experimental time scale (Fig. [3]^b)), 
in perfect agreement with the usual assertion that the 
system becomes glassy around = 0.58 [3J- We also ob- 
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FIG. 4. Probability distribution of neighbor changing events. 
The number of events decreases exponentially with the num- 
ber of changed neighbors in the event. The solid curve is an 
exponential fit. Inset: The number of events over a given 
time period varies as a power law with the free energy change 
associated with the events. 

serve rather large fluctuations in the time dependence of 
free energy suggesting that every now and then the sys- 
tem becomes trapped in local free energy minima (Fig. 
[3](d)) , as is also expected for a glass [2711125] . 
To look at the aging in more detail, we investigate 
small subsystems, so that the fluctuations are more pro- 
nounced. We calculate the free energy for many small 
boxes, each containing around 60 particles, in the col- 
loidal glass at cj) = 0.60 (Fig. [3^c)). The free energy of 
the glass is found to be very heterogeneous: The free en- 
ergy varies strongly with the position of the box, e.g. the 
free energy may change about 20%, or more than lfcgT 
per particle from one region to another at any given time. 

We are now in the position to address the central ques- 
tion about glassy dynamics, namely how the aging of the 
free energy is related to the rearrangements of particles. 
To investigate this, we computed the number of changed 
nearest neighbors for each pairs, located in regions that 
contain around 200 particles, from the Voronoi decompo- 
sition [52]. Comparing the neighbors of all the particles 
after each time step allows us to detect rearrangements 
as events for which a certain number of neighbors change 
between two time steps (Fig. [4]). We find that the prob- 
ability of occurrence of a neighbor changing event de- 
creases exponentially with the magnitude of the event, 
i.e. the number of changed neighbors. If we take AF to 
be the absolute value of the free energy difference between 
two consecutive snapshots (see Fig. (3^d)) averaged over 
events with the same number of changed neighbor pairs, 
we find that the number of events obeys a power law 
scaling in terms of the free energy changes with a power 
~ —1.2 ± 0.15 (Fig. [I]). Surprisingly this is analogous to 
the Gutenberg-Richter law [30] in seismology where the 
probability for large earthquakes to happen decreases ex- 
ponentially with the magnitude of the earthquakes, the 
latter being logarithmic in terms of the released energy; 
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in this case number of events is a power law in terms of 
the released energy with a power ~ — 1 ± 0.4 [31]. 
In summary, we present the first measurements of the free 
energy of a glassy system. The concept of the free energy 
landscape and its deep minima is central in our current 
theoretical understanding of glassy systems. However, a 
direct measurement was so far not possible but has been 
achieved here and reveals the strong heterogeneity of the 
free energy and its link to the slow structural relaxation 
that is the hallmark of glassy behavior. 
We would like to thank the Stichting voor Fundamenteel 
Onderzoek der Materie (FOM) and Shell for the financial 
support of the present research work. 
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